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Abstract 

For pairs u, p, of density operators on a finite dimensional Hilbert 
space of dimension d I call /c-fidelity the d — k smallest eigenvalues 
of \i/ujy/p\. A;- fidelities are jointly concave in uj,p. This follows from 
representing them as infima over linear functions. For k = known 
properties of fidelity and transition probability are reproduced. Partial 
fidelities characterize equivalence classes of pairs of density operators 
which are partially ordered in a natural way. 

1 Introduction 

For two unit vectors, ip and (p, of an Hilbert space the quantity \(ip,ip)\ 2 
is their transition probability. It is the squared modulus of their transition 
amplitude, (tp, <p). Assume the state of the quantum system is \if))(ip\. A von 
Neumann measurement, designed to decide whether the quantum system is 
in the state \(p)((p\, prepares this state with probability tp)\ 2 - 
Notice further that two pairs of unit vectors are unitarily equivalent iff they 
enjoy equal transition probabilities. 

All that becomes more complex, ||, if two density operators, pi and p 2 , 
are considered on a Hilbert space TC, and the quantum system is in state, say, 
p\. The algebra of operators on TC will be called £>. One can choose vectors 
ipj in the direct product TC <S> TC such that 

TrApj — {ipj, (A <g> A e B, j = 1,2 (1) 

The transition probability between ifii and ^2 is not determined by the pair 
Pi, pi. But running through all possible arrangements ([l|), the numbers 
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|(V'2)'0i)| 2 fill completely an interval [0,p] of real numbers. The largest one, 
the upper bound of this interval, is called transition probability between pi 
and P2 and is denoted by P(pi,p2). Thus a von Neumann measurement 
in H Cg> 7i can cause a transition p\ \— > p 2 with a probability bounded by 
P{px,P2)- The bound can be reached by suitable measurements in the larger 
system. 

Now I call attention to possibilities to characterize P intrinsically, i. e. with- 
out leaving the quantum system in question. The first one comes rather 
directly from (|l|). Let us call transition functional from p\ to p2 every linear 
functional on B of the form 

A — >Tr uA:= (2) 

which arises from a setting ([[]). The operators v may be called transition 
operators from p\ to p 2 . Generally, v is not Hermitian: Exchanging the roles 
of and i\)2 the operator v becomes its Hermitian adjoint, v* . 
Now (0) is a transition functional if and only if 

| Tr AivAl | 2 < {TrA lPl A* 1 ) (Tr A 2 p 2 A* 2 ), AeB (3) 

and it follows from the definition of P that 

P(P2, Pi) = m & x I Tr v | 2 (4) 

where one takes the maximum over all transition operators from p\ to p 2 - 
Calculating the maximum in (§) is a standard exercise with a well known 
outcome. Before writing it down I would like to explain the following: 
The transition probability is separately concave in every one of its arguments. 
However, taking the root of P, the concavity properties become dramatically 
enhanced: \fP is jointly concave 0. In the following the square root of 
the transition probability will be called fidelity and will be denoted by F, 
essentially following a proposal of Richard JozsaQ. Thus 

F(u, p) := y/P(u,p) = Tr^/W 72 ) 172 (5) 
The assertion that F is jointly concave is seen from 

F(u,p) = ^inf(Tr Au + TiA^p), A > 0, A' 1 6 B (6) 

1 Jozsa introduced the word fidelity for the transition probability. Its present usage is 
not unique. I think the peculiar properties of \/~P need an extra notation anyway. 
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which is the finite dimensional version of a representation of \fP = F as 
an infimum of linear functionals, valid for pairs of states on von Neumann 
and on C*-algebras, see ||. The representation is related to another one of 
equal generality estimating P(u, p) from above by the product of TruM and 
TrpA^ 1 , with A an invertible positive operator, see @ for a partial result 
and || for the C*-case in full generality For finite dimensions these well 
know results are reproduced by setting k = in the expressions fll3|) and 
(0) below. 

As a matter of fact, the equality of F (or of P) for two pairs of density 
operators do not imply their unitary equivalence. This pleasant feature, valid 
for pure states, is missing for mixed ones. Looking at one may wonder 
whether it is not possible to get a whole series of concave invariants by taking 
other suitable sets of operators than the invertible positive operators in the 
expression (||). To give an affirmative answer belongs to the issues of the 
present paper. By the partial fidelities one gets a reasonable classification of 
pairs of density operators, coarser than unitary equivalence would give. 

All what follows remain in finite dimensions. By modifying certain set- 
tings and by adding new arguments, Peter M. Alberti, 0, was able to extend 
essential parts of what follows to von Neumann Algebras. His results are par- 
ticularly satisfying for type IIx- 

2 /c-Fidelities 

Let H be a finite-dimensional Hilbert space and d = dim 7-^. 
The spectrum, spec(^4), of an operator A is the family of roots of the polyno- 
mial det(A — Al) counted with their correct multiplicities. If the spectrum is 
real we assume the set spec(yl) decreasingly ordered. This convention applies 
to every diagonalizable operator with real eigenvalues and in particular to 
every Hermitian one. Consider now 

spec(( v / ^Pv / ^) 1/2 ) = spec(( v ^a; v ^) 1 / 2 ) = { X 1 > A 2 >...,> X d } (7) 

so that, according to @, the sum of the lambdas is the fidelity. The spectrum 
(0) is equal to the ordered singular numbers of \fu^[p and of ^Jp~\fuj. 
I define partial fidelities simply by summing up parts of the spectrum (0). 
Definition: For < k < d — 1 

F k (u;,p):=J2^, fc = 0,l,...,d-l. (8) 

j>k 
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If k > d then F k = 0. For the time being F k will be called k-th partial fidelity, 
or simply k -fidelity of the pair uo and p. 

An important point is to add: I do not necessarily require that p and uj 
have trace one. Indeed, on a finite dimensional Hilbert space (§) is naturally 
defined for pairs from the cone of positive operators and 

y/cF k {uj,p) = F k (cuj,p) = F k (uj,cp), c>0 (9) 

for positive real numbers c. Notice the properties: 

a) F k is symmetric in its arguments. 

b) F is just the fidelity F. 

c) For pairs of pure density operators it is F k = for k > 0. 

d) If F k ^ then rank(kj) > k and rank(p) > k necessarily. 

e) F k is unitarily invariant, i.e. invariant by the simultaneous transformation 
p -> UpU* , uj -> t/wC/*. 

However, a deeper justification for the definition above is in 
Theorem 1 

The partial fidelities are concave functions of the pairs {u, p} : 

J2Pj F k(^j,Pj) < F k (J2PjUj,J2PiPi) ( 10 ) 

3 3 i 

for any probability vector pi,p 2 , . . . and arbitrary pairs {ui, p^. □ 
The theorem is a consequence of a new relation representing F k as an infimum 
of linear functionals quite similar to (|6]). It estimates partial fidelities linearly 
as close as possible from above. To get the announced representation I am 
going to define the set PAIRS which consists of all pairs {A, B} of positive 
Hermitian operators, A, B, such that 

ABA = A, BAB = B (11) 

Let {A, B} be such a pair. It follows (AB) 2 = AB immediately. Because 
Q = AB is a product of two positive operators it is diagonalizable. On the 
other hand we see Q 2 = Q so that its spectrum consists of zeros and ones. 



Therefore, the trace of Q is equal to the rank of Q. Now (11) says QA = A 
and BQ = B implying that the ranks of A and B cannot be larger than the 
rank of Q. Now Q = AB shows that neither the rank of A nor the rank of 
B can be less than that of Q. Altogether we have: 
Lemma 1: For all {A, B} E PAIRS 

rank(y4) = rank(S) = mnk(AB) = Tr AB (12) 
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is an integer called rank of the pair {A, B}. □ 

Definition: PAIRS m consists of all pairs from PAIRS of rank m. 

The promised representation of the fc-th fidelities is in the following theorem. 

Theorem 2 

Let m + k = dim7i. Then 

F k (u,p) = -ird(TrAou + TrBp), {A, B} G PAIRS m . (13) 

One can deduce from (|T3|) the following inequality: 

F k (u,p) 2 = inf(Tr Au^TrBp), {A, B} G PAIRS m . (14) 

The point is that with {^4, B} also {XA, X~ 1 B} is contained in PAIRS m for 
A > 0. After this trivial substitution, the right hand side o f (O ) is of the 
form Xa + X~ 1 b. Taking the infimum of over A results in 2yab, and ( |I~4"D is 
derived from ([13]). (|14]) , suitably reformulated, is known on C*-algebras if 
k = 0, see §. 

Theorem 1 is a consequence of theorem 2. We shall prove theorem 2 in 
the next section, at first assuming p invertible, (which would be sufficient 
for theorem 1). Then, by continuity arguments, we can allow for all u and 
p. However, before going into the proof, we have to look at a "hidden" 
symmetry of the ^-fidelities. 

3 The symmetry group of the /c-fidelities 

Let us denote by T the multiplicative group of all invertible operators acting 
on 7i. With X G T we define the X -transform of a pair {u, p} by 

{u,p} x -^{XuX^iX-'ypX- 1 } (15) 

The transformations create orbits of T in the set of pairs. Two pairs, {uj, p} 
and {u)',p'}, are called V -equivalent iff there is X G V such that {a/, p'} is 
the X-transform of {u, p}. 

Lemma 2: The fc-fidelities of T-equivalent pairs are equal for every k 

F k (u, p) = F k (u/, p') if W, p'} = {u, p} x (16) 
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For the proof we start with the identity 

(XluX*)(X- 1 )* pX- 1 = XupX- 1 

saying that the spectrum of up is an invariant for T-equivalent pairs. It 
suffices to show that 

spec(^p) = spec(pw) = { A?, \\, . . . , A^ } (17) 

follows from (0). By ^/p(up)( y /p)~ 1 = ^fpoj^fp this is true for invertible p. 
The assumption of invertibility can be removed by continuity. □ 
Substituting in (pD p = X*rX we may rewrite (|16|) as 

F k (u, X*rX) = F k (XuX*, t) (18) 

In ([TSj) , by continuity, also X need not be invertible. 

Only for the purpose of the the following proof we abbreviate the right- 
hand-side of fll3T) by G k (u,p), 

G k (u,p) := - mf^TiAuj + TrBp), {A, B} G PAIRS m . 

We observe that also G k allows for T-invariance. Indeed, with a pair {^4, B} 
its transform 

{A, B} x := {X*AX, X- l B(X- 1 )*} (19) 
is also in PAIRS m , and we get the trace identities 



whenever 
Therefore, 



TruM = Ttu'A', TrpB = Tip'B' 
{A',B'} = {A,B} X , W,p?} = {u,p} x 



G k (uj,p) = G h {u)',p') if {uj',p} = {u,p} x (a) 

Now we start the proof of theorem 2. As we just have seen, both sides of (p~3|) 
do not change along a T-orbit of X-transforms (JT5D- 

Step 1 in the proof is to show F k < G k . To this end we first assume a pair 
{oo,p} with an invertible p. We transform the given pair according to (15) 
by X = y/p~. The new pair is {uj 1 ', 1} with u' = ^fpuj-^fp. By (a) and lemma 
2 it suffices to prove F k < G k for the new pair, i. e. we estimate G k (uj', 1) 
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from below. 

We choose a pair {A, B} from PAIRS m arbitrarily. Let fa,fa,... be an 
eigenbasis of A and a\, a-i, ■ ■ ■ the corresponding eigenvalues. By sandwiching 
ABA = A between these eigenvectors of A one gets 

(fa, Bfa) = a^dik for i,k>m 

Now we can write 

m m 

Tr Aw' + Tr B = £ a^fa u'fa) + £ af 1 + £ (</>„ 50,) 

1 1 j>m 

With positive reals a and x it holds ax + a -1 > 2-y/ir. Using that inequality 
to estimate the first two sums from below and neglecting the last term, we 
arrive at 

m 

TtAlo' + TtB > 2Y,\/(fa,u'fa) 
i 

The square root is concave. Hence, see [fTCA , equ. 1-46, 

m m 

£ \l (fa "'fa) > J2(fa, Vu/fa) > F k (u, p) 
1 1 

The last inequality sign is an estimation of the m smallest eigenvalues (due 
to Fan and Horn) and respecting F k (uj' , 1) = F k {oJ, p). It results 

- inf(TrAa; + Tr5p) > F k (u),p) (20) 

at first for pairs {uj' , 1} and then, by T-invariance, for all pairs {uj, p} with 
invertible p. However, both sides of (pOp are continuous in u and p. Thus 
step one terminates in the validity of fl2~0|) for all uj and all p. The inequality 
is equivalent to G/. > F k . 

In step 2 we show G k < F k at first for invertible p. As above we reduce 
the problem by T-invariance to that of a pair consisting of u' and 1. We 
now choose fa, . . . , (j) m to be eigenvectors of uj' belonging to the m smallest 
eigenvectors of uj' . The latter are A| +1 , . . . , A| +m . by lemma 2. Define 

m m 

i i 
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and consider 



m m 



TtA'uj' + TrB> = £ a 3 X) +k + £ a- 



-i 



j=l 1 



If Xj + k > we choose aj = Xj+ k . Otherwise we set a 3 —c 1 > arbitrarily. 
If n of the m eigenvalues Xj+k are zero, then our convention implies 



and, hence, Gk < Fk + nc. Since c can be made arbitrarily small we arrive at 
the wanted inequality Gk < Fk. Now, relying on T-invariance ([Tj]) and (a), 
the inequality is shown true for all pairs of invertible density operators. 
Combining step one and two we see: F k (uj,p) = Gk(uJ,p) if both arguments 
are invertible. Hence F% is concave for these pairs. But Fk is a continuous 
function of uj and p by (§). Therefore, Fk is jointly concave and theorem 1 
is valid. 

But one knows that a concave function is semi- continuous from below, see 
J7J, theorem 10.2, where semi-continuity from above is stated for convex 
functions. Because Fk is continuous and concave it dominates every function 
which is concave and coincides for convexly inner points with Fk- This means 
Fk > Gk always. Now step one of the proof provides Fk = Gk- □ 

4 Equivalence and partial order 

It is tempting to collect pairs of positive (density) operators into equivalence 
classes according to their partial fidelities. For the purpose of the present 
paper we call two pairs equivalent, and we write 



iff their /c-fidelities are equal, Fk(u,p) = Fk(u',p') for k — 0, 1, . . . , d — 1. 
The relation ~ is an equivalence relation. Notice that {uj,p} ~ {p,uj,}. 
Generally, an equivalence class contains a lot of T-orbits. But there is an 
important exception: 

Lemma 3: If both operators, uj and p, are invertible, the equivalence class 
of {uj, p} consists exactly of all pairs {uj, p} x , X £ V .□ 

Proof: The assumption is valid if and only if does not belong to the 
eigenvalues (0). This takes place if the smallest one is different from zero, 



rn 



TtA'uj' + TtB' = 2Y, ^j+k + nc = 2F k {uj 1 , 1) + nc 



{u,p} ~ {uj',p} 
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hence iff F^-i 7^ 0. Hence, if the assumption of the lemma is valid for one 
member of an equivalence class, then it is true for all members. Let {u\,p\\ 
be in the equivalence class of {u>,p}. Transforming the latter by X = ^fp 
and the former by X\ = ^fp~[ by the receipt ([0J) results in accordingly T- 
equivalent pairs {uj', 1} and {ou[, 1}. Being in the same equivalence class, u' 
and oj[ have to have equal eigenvalues and they are even unitarily equivalent. 
Thus all the pairs considered belong to the same T-orbit. □ 

Let us write {u',p'} < {uj, p} if both, u> — co 1 and p — p', are positive 
operators. A simple example is as follows: Write uj = uj' + uj , p = p' + p , 
and assume orthogonality between uj' and p and between p' and u , i. e. 
u>op' = 0, pqoj' = 0. Then {uj,p} and {u',p'} belong to the same equivalence 
class. To see what we can learn from {uj',p'} < {uj,p} generally, we proceed 
in two steps, {u/, p'} < {uj,p'} and {uj,p'} < {u,p}. Consider the second 
one. It implies ^Jujp'^/uj < ^Jujp^To and, because taking the square root does 
not destroy the inequality, 

(v^p'v^) 1/2 < (v^ Pv ^) 1/2 

The sums of its m smallest eigenvalues, which are the partial fidelities, obey 
the same inequality. Further, if the traces of both positive operators happen 
to be equal, the operators themselves have to be equal one to another. In 
repeating the arguments for the first step and combining both, we arrive at 
Lemma 4: If 

{u,p}>W,f/} (22) 

then 

F k (u,p)>F k (u',(/), fc = 0,l,...,d-l (23) 

If in addition to ( [22] ) F(u>,p) = F(u',p') is true, then all partial fidelities 
must be equal in pairs, and the two pairs belong to the same equivalence 
class: {u, p} ~ {uj', p'}. □ 

Given uj, p, Alberti || has shown, even in the C*-category, that there is 
one and only one pair {uq,Pq} which enjoys the same transition probability, 
(and, therefore, the same fidelity), and which is minimal with respect to >. 
This minimal pair satisfies 

{^o,Po} < W,p'} 

whenever 

W,p'} < W,p} and F(u',p) = F(u,p) 
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is valid. 

We see that every equivalence class contains a minimal pair and, therefore, a 
T-orbit of minimal pairs. It is tempting to believe that there is only one min- 
imal T-orbit in every equivalence class of pairs. But I do not know whether 
that conjecture is true. 

Now one may go a step further, anticipating the ideas of majorization, 
||, or those of partially ordering orbits belonging to certain classes of trans- 
formations, [fLOfl . To do so, let us call {coi, p{\ F-dominated by {u>2, P2} iff 



F k {u uPl )<F k {u 2 ,p 2 ), A; = 0,1,2,... (24) 
From theorem 1 we get the following 

Corollary: If {102, P2} is contained in the convex hull of the ~equivalence 
class of {ui,pi} then ( |24"D takes place. □ 

We thus get a new partial ordering (or majorization tool) for pairs of 
positive (density) operators which seems worthwhile to investigate. There 
is a link, indeed a morphism, to singular number majorization. Denote by 
sing(.B) the decreasingly ordered singular numbers of the operator B, that is 

sing(S) = spec("\/ B*B) = spec(V BB*) = sing(S*), 

and by sing[S] the set of all operators C with sing(C) = sing(S), the singular 
number class of B. In particular 

spec ((y/upy/u) 1/2 ) = sing(y'pV^) 

There are many useful rules governing the partial order of the singular num- 
ber classes, see ||, 9.E and [10], 2.4 (theorem 2-8), 2.5. With them one easily 



proves 

Lemma 5: The following items are mutually equivalent: 

a) {ujx, pi} is F-dominated by {^2,^2}- 

b) Tj^yffh is contained in the convex hull of singf^/aJi^/pi]. 

c) There are finitely many operators Aj, Bi, all with operator norms not 
exceeding 1, such that 

UJ2\/p~2 = Ay/uhy/plBi 
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5 More about PAIRS 



It is our aim to get some insight into the structure of PAIRS. Let {A, B} e 
PAIRS m with < m < d = dim Ti. k is defined by k + m = d. Let us write 
A, B as block matrices with respect to an eigenvector basis of A as in the 
proof of theorem 2. Then, with is a positive m x m matrix An, 

An 0\ „ / Bn B12 



Here Bn is m x m, B i2 is m x k, and B 22 is k x k. The equation ABA = A 
results in Bn = An - Having this in mind, one gets from BAB = B 

Bn = A u x , B 2 2 = B 2 iAiiBi 2 , B\ 2 = B 2 i (26) 

Notice that B u can be chosen arbitrarily: Given the first member, A, of the 
pair, B depends freely on km complex parameters. 

There is a further representation of the pairs in PAIRS m . Call 2m vectors, 
■01, ij>2, ■ ■ ■ , ipm, <Pi, <fi2, ■ ■ ■ , fm, a bi-orthogonal system of length m if 

(i/ji,(Pj) = Sij, i,j = 1,2,..., m (A) 

Together with m positive numbers, ai, a 2 , ■ ■ ■ , a m , we obtain from (A) a pair 
of operators 

A = J2a l \^)(ip i \, B = Y< a 7 1 Wi){<Pi\ (B) 

for which ABA = A and BAB = B can be checked. Let us prove that every 
pair from PAIRS m can be gained by this procedure. 

Let {A, B} G PAIRS m . Because AB is diagonalizable with eigenvalues 
and 1, there is X e T such that XABX^ 1 is a Hermitian projection operator. 
Hence the operators 

Ai = XAX*, Bi = (X^YBX- 1 
commute. Therefore there is a representation 

A 1 = ^2a i \<f> i )((/> i \, Bi = J2 a i 1 \ ( Pi)( ( t ) i\ 
with m orthonormal vectors <pi, . . . , <p m . But 
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is bi-orthogonal with length m. Transforming A\ and Bi back to A and B 

gives the desired representation of the pair. 

The bi-orthonormal system (A) of (B) can be chosen balanced : 

(ipi,tpi) = ((puipi), i = l,...,m (C) 

Indeed, the necessary changes in the norms can be compensated by adjusting 
the aj. Now we insert (B) into the right hand side of (|13|) and observe 

m 

Tr Ato + Tr Bp = ^(oj < ij)j\u)\ipj > < <fj\Q\fj >) 

1 

By varying the free parameters dj we arrive at 
Theorem 3 

Let m + k = dimH. Then 



F k (u,p) =M'*r f J(il>i,uj'il) i )((pi,p(pi) (27) 



i=i 



where the infimum runs through all balanced bi-orthogonal systems of length 
m. □ 

Finally, assume the infimum in (|T3|) is attained by {A, B} 6 PAIRS m , 

F k (u,p) = ^(Tr^w + TrSp), m + k = dimH (28) 
If we vary the minimizing pair, the first variation must vanish, 

(;r) A s ) » c ( s ) = ( Tr A * u + Tr b *p) 
where, with X s = exp sY and any operator Y, 

A s = X* S AX S , B s = X^BiX*)- 1 
We perform the first derivative and obtain 

(pj A s = Y* A + AY, (pj B s = —YB — BY* 

After inserting in c(0) = and an rearrangement it results 
TrY(Ap-ojB)* + TrY*(Ap - ujB) = 

12 



As Y could be chosen arbitrarily, we arrive at 



Ap = uB (29) 

as a necessary condition for the validity of (2~S). 

Is there {A, B} 6 PAIRS m fulfilling (^) and minimizing (|T3|)? If we can 
r-transform u>, p, to the form r, r, we are done. Indeed, we then can choose 
a projection operator P m onto the m smallest eigenvalues of r and we get 

F k (r, t) = Tr P m r, {P m , P m } G PAIRS m 

i. e. the problem is solved in that case. Now, if uj and p are both invertible, 
there is a unique positive X such that 

XuX = X^pX- 1 :=r, X > (30) 

The choice ([KJ) ensures @ with A = XP m X, 5 = r 1 ?^ 1 . To get X 
one has to solve X 2 uX 2 = p. There is a unique positive solution X which is 
the square root of the geometric mean between p and u^ 1 . 

X 2 = UJ- 1 /\U 1 / 2 P UJ 1 I 2 ) 1I2 LJ- 1 / 2 

as one can convince oneself by inserting into X 2 ujX 2 = p. 
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